clear all
load W_time1_age0_eh1_comp11.out
n = 50;
w_f_g=W_time1_age0_eh1_comp11; % Read wages, g and f distributions
wf_1=w_f_g(:,1);        % Wage in formal sector
wi_1=w_f_g(:,2);        % Wage in informal sector
wf_2=w_f_g(:,7);      % Wage in formal sector
wi_2=w_f_g(:,8);      % Wage in informal sector

load Rw_1_f_gamma_0_GE Rw_1_f_gamma_0_GE
load Rw_2_f_gamma_0_GE Rw_2_f_gamma_0_GE
load Rw_1_i_gamma_0_GE Rw_1_i_gamma_0_GE
load Rw_2_i_gamma_0_GE Rw_2_i_gamma_0_GE

% Individual model (men)
load Rw_1_f_ind
load Rw_1_i_ind

x_head = wi_1; x_spouse = wi_2;

plot(x_spouse,Rw_1_i_gamma_0_GE,'-b',x_spouse,Rw_1_f_gamma_0_GE,'-g',x_spouse,Rw_1_i_ind*ones(50,1),'--b',x_spouse,Rw_1_f_ind*ones(50,1),'--g','LineWidth',2.5)
 
xlabel('Wage of women','FontSize',15)
ylabel('Reservation wage of men','FontSize',15)
box off
 
legend('$\hat{w}_{ni->ii}$','$\hat{w}_{ni->fi}$','$\hat{w}_{n->i}$','$\hat{w}_{n->f}$','Interpreter','latex','FontSize',20)


print -dpdf Rw_ind_hh_search_men.pdf



